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Abstract — hypergraph  is  a  set  V  of  vertices  and  a  set  of 
nonempty  subsets  of  V,  called  hyperedges.  Unlike  graphs,  hy¬ 
pergraphs  can  capture  higher-order  interactions  in  social  and 
communication  networks  that  go  beyond  a  simple  union  of  pair¬ 
wise  relationships.  In  this  paper,  we  consider  the  shortest  path 
problem  in  hypergraphs.  We  develop  two  algorithms  for  finding 
and  maintaining  the  shortest  hyperpaths  in  a  dynamic  network 
with  both  weight  and  topological  changes.  These  two  algorithms 
are  the  first  to  address  the  fully  dynamic  shortest  path  problem 
in  a  general  hypergraph.  They  complement  each  other  by  parti¬ 
tioning  the  application  space  based  on  the  nature  of  the  change 
dynamics  and  the  type  of  the  hypergraph.  We  analyze  the  time 
complexity  of  the  proposed  algorithms  and  perform  simulation 
experiments  for  random  geometric  hypergraphs,  energy  efficient 
routing  in  multichannel  multiradio  networks,  and  the  Enron  email 
data  set.  The  experiment  with  the  Enron  email  data  set  illustrates 
the  application  of  the  proposed  algorithms  in  social  networks 
for  identifying  the  most  important  actor  and  the  latent  social 
relationship  based  on  the  closeness  centrality  metric. 

Index  Terms — Dynamic,  hypergraph,  hyperpath,  shortest  path. 


1.  Introduction 

A  GRAPH  is  a  mathematical  abstraction  for  modeling  net¬ 
works,  in  which  nodes  are  represented  by  vertices  and 
pairwise  relationships  by  edges  between  vertices.  A  graph  is 
thus  given  by  a  vertex  set  V  and  an  edge  set  E  consisting  of 
cardinality-2  subsets  of  U .  A  hypergraph  is  a  natural  extension 
of  a  graph  obtained  by  removing  the  constraint  on  the  cardinality 
of  an  edge:  Any  nonempty  subset  of  V  can  be  an  element  (a  hy¬ 
peredge)  of  the  edge  set  E  (see  Fig.  1).  It  thus  captures  group 
behaviors  and  higher-dimensional  relationships  in  complex  net¬ 
works  that  are  more  than  a  simple  union  of  pairwise  relation¬ 
ships.  Examples  include  communities  and  collaboration  teams 
in  social  networks,  document  clusters  in  information  networks, 
and  cliques,  neighborhoods,  and  multicast  groups  in  communi¬ 
cation  networks. 
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Fig.  1.  Example  hypergraph  with  four  hyperedges:  (vi,  V2,  V3,  vq), 

(u2,  U3,  U4,  U5),  (U6,  U7,  V8,Vq),  and  (ub,  Ug,  Ug). 


While  the  eoneept  of  hypergraph  has  been  around  sinee  the 
1920s  (see,  for  example,  [1]),  many  well-solved  algorithmie 
problems  in  graph  remain  largely  open  under  this  more  general 
model.  Here,  we  address  the  shortest  path  problem. 

A.  Shortest  Path  Problem  in  Graphs 

The  shortest  path  problem  is  perhaps  one  of  the  most  basie 
problems  in  graph  theory.  It  asks  for  the  shortest  path  between 
two  vertiees  or  from  a  souree  vertex  to  all  the  other  vertiees  (i.e., 
the  single-souree  version  or  the  shortest  path  tree).  Depending 
on  whether  the  edge  weights  ean  be  negative,  the  problem  ean  be 
solved  via  Dijkstra’s  algorithm  or  the  Bellman-Ford  algorithm. 

The  dynamie  version  of  the  shortest  path  problem  is  to  main¬ 
tain  the  shortest  path  tree  without  reeomputing  from  serateh 
during  a  sequenee  of  ehanges  to  the  graph.  A  typieal  ehange  to  a 
graph  ineludes  weight  inerease,  weight  deerease,  edge  insertion, 
and  edge  deletion.  The  last  two  types  of  ehanges  model  network 
topologieal  ehanges,  but  they  ean  be  eoneeptually  eonsidered 
as  speeial  eases  of  weight  ehanges  by  allowing  infinite  edge 
weights.  Thus,  if  the  sequenee  of  ehanges  eontains  only  edge 
deletion  and  weight  inerease,  we  eall  it  a  deeremental  problem; 
if  it  eontains  only  edge  insertion  and  weight  deerease,  we  eall 
it  an  ineremental  problem.  Otherwise,  we  have  a  fully  dynamie 
problem.  If  multiple  edges  ean  ehange  simultaneously,  then  it  is 
a  bateh  problem.  Example  dynamie  shortest  path  algorithms  for 
graphs  ean  be  found  in  [2]-[5]. 

B.  Shortest  Path  Problem  in  Hypergraphs 

Both  the  statie  and  dynamie  shortest  path  problems  have  a 
eorresponding  version  in  hypergraphs.  A  hyperpath  in  a  hyper¬ 
graph  is  a  sequenee  of  hyperedges  with  two  adjaeent  hyperedges 
sharing  at  least  one  vertiees.  The  weight  of  a  hyperpath  is  the 
sum  of  the  weights  of  eaeh  hyperedge  on  this  hyperpath.  The 
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shortest  hyperpath  problem  in  weighted  hypergraphs  where  a 
nonnegative  weight  is  attaehed  to  eaeh  hyperedge  is  thus  elearly 
defined.  The  statie  shortest  hyperpath  problem  was  eonsidered 
by  Gallo  et  al.  [6]  in  whieh  Dijkstra’s  algorithm  was  extended 
to  obtain  the  shortest  hyperpaths.  Ausiello  et  al.  proposed  a 
semi-dynamie  minimum  length  hyperpath  algorithm  with  time 
eomplexity  0{K\H\  log  h-\-Q)  for  direeted hypergraphs.  While 
they  eonsidered  a  more  general  definition  of  hyperpath  length 
(whieh  ean  be  an  arbitrary  funetion  of  the  hyperedge  weight  on 
the  hyperpath),  their  dynamie  algorithm  was  developed  under 
a  very  restrietive  assumption.  Speeifieally,  it  was  assumed  that 
the  length  of  a  path  ean  only  eome  from  a  finite  diserete  set  with 
size  K.  Sinee  the  time  eomplexity  of  their  algorithm  is  linear  in 
K,  their  algorithm  is  not  applieable  to  hypergraphs  where  edge 
weights  ean  be  any  real  number  as  eonsidered  in  this  paper.  Fur¬ 
thermore,  they  eonsidered  only  the  ineremental  part  of  the  dy¬ 
namie  maintenanee,  while  in  this  paper  we  eonsider  both  inere¬ 
mental  and  deeremental  problems.  A  dynamie  algorithm  for  the 
bateh  problem  for  a  speeial  elass  of  hypergraphs  was  developed 
in  [3], 

With  the  exeeption  of  the  above  few  studies,  the  shortest  hy¬ 
perpath  problem  remains  largely  unexplored.  To  the  best  of  our 
knowledge,  no  algorithms  exist  for  the  fully  dynamie  problem 
in  a  general  hypergraph. 

In  this  paper,  we  develop  two  fully  dynamie  shortest  path  al¬ 
gorithms  for  general  hypergraphs.  These  two  algorithms  eom- 
plement  eaeh  other,  with  eaeh  preferred  in  different  types  of  hy¬ 
pergraphs  and  ehange  dynamies. 

Referred  to  as  the  HyperEdge-based  Dynamie  Shortest  Path 
algorithm  (HE-DSP),  the  first  algorithm  is  an  extension  of  the 
dynamie  Dijkstra’s  algorithm  for  graphs  to  hypergraphs  (par¬ 
allel  to  Gallo’s  extension  of  the  statie  Dijkstra’s  algorithm  to 
hypergraphs  in  [6]).  The  extension  of  the  dynamie  Dijkstra’s  al¬ 
gorithm  to  hypergraphs  is  more  involved  than  that  of  the  statie 
Dijkstra’s  algorithm.  This  is  due  to  the  loss  of  the  tree  strueture 
(in  the  original  graph  sense)  in  the  eolleetion  of  the  shortest  hy¬ 
perpaths  from  a  souree  to  all  other  vertiees.  Sinee  the  dynamie 
Dijkstra’s  algorithm  relies  on  the  tree  strueture  to  update  the 
shortest  paths  after  an  ineremental  ehange,  speeial  eare  must  be 
taken  when  extending  it  to  hypergraphs. 

The  seeond  algorithm  is  rooted  in  the  idea  of  Dimension  Re- 
duetion  and  is  referred  to  as  DR-DSP.  The  basie  idea  is  to  re- 
duee  the  problem  to  that  in  an  indueed  graph  derived  from  the 
original  hypergraph.  The  indueed  graph  of  a  hypergraph  has  the 
same  vertex  set  and  has  an  edge  between  two  vertiees  if  and  only 
if  there  is  at  least  one  hyperedge  eontaining  these  two  vertiees 
in  the  original  hypergraph.  The  weight  of  an  edge  in  the  indueed 
graph  is  defined  as  the  minimum  weight  among  all  hyperedges 
eontaining  the  two  vertiees  of  this  edge.  The  shortest  hyperpath 
in  the  hypergraph  ean  thus  be  obtained  from  the  shortest  path  in 
the  indueed  graph  by  substituting  eaeh  edge  along  the  shortest 
path  with  the  hyperedge  that  lent  its  weight  to  this  edge.  The 
eorreetness  and  advantage  of  this  algorithm  are  readily  seen: 
The  definition  of  weight  in  the  indueed  graph  eaptures  the  min¬ 
imum  eost  offered  by  all  hyperedges  in  ehoosing  a  path  between 
two  vertiees,  thus  ensuring  the  eorreetness  of  the  algorithm;  the 
reduetion  of  a  hypergraph  to  its  indueed  graph  removes  many 
hyperedges  from  eonsideration  when  finding  the  shortest  path, 
leading  to  effieieney  and  agility  to  dynamie  ehanges. 


As  will  be  shown  in  the  time  eomplexity  analysis  given  in 
Seetion  V,  HE-DSP  is  more  efiieient  in  hypergraphs  that  are 
densely  eonneeted  through  high-dimensional  hyperedges  and 
for  network  dynamies  where  ehanges  often  oeeur  to  hyper¬ 
edges  that  are  not  on  the  eurrent  shortest  hyperpaths.  DR-DSP 
has  lower  eomplexity  when  hyperedge  ehanges  often  lead  to 
ehanges  in  the  shortest  hyperpaths.  This  is  usually  the  ease  in 
networks  where  hyperedges  in  the  shortest  hyperpaths  are  more 
prone  to  ehanges  due  to  attaeks,  frequent  use,  or  higher  priority 
in  maintenanee  and  upgrade.  While  we  foeus  on  undireeted 
hypergraphs  in  the  paper,  the  two  algorithms  apply  to  direeted 
hypergraphs  as  diseussed  in  Seetion  VI.  They  also  apply  to 
bateh  hyperedge  ehange  problems  with  minor  modifieations. 

C.  Applications 

Shortest  path  eomputations  on  hypergraphs  ean  be  applied  to 
eommunieation  as  well  as  soeial  networks.  An  example  appliea- 
tion  is  routing  in  multiehannel  multiradio  ad  hoe  networks  [8]. 
As  detailed  in  Seetion  VII,  we  eonsider  energy-effieient  routing 
where  the  link  metrie  is  energy  eonsumption  ineluding  both 
transmitting  and  reeeiving  energy.  To  eapture  the  presenee 
of  different  neighbor  sets  assoeiated  with  different  ehannels 
available  to  a  partieular  node,  a  direeted  hypergraph  model 
is  needed  in  whieh  a  hyperedge  is  formed  from  a  vertex  to 
its  neighbors  that  shares  a  partieular  ehannel  and  different 
ehannels  available  to  this  vertex  lead  to  different  hyperedges. 
The  weight  of  a  hyperedge  is  given  by  the  total  energy  eost 
of  a  transmission  from  the  sending  vertex  using  the  ehannel 
eorresponding  to  this  hyperedge.  A  shortest  hyperpath  from  the 
souree  to  the  destination  is  the  least-energy-eost  route,  and  as 
network  topology  ehanges,  a  dynamie  algorithm  is  required  to 
maintain  the  shortest  hyperpath. 

In  soeial  networks,  information  (results,  event  reports,  opin¬ 
ions,  rumors,  ete.)  propagates  through  diverse  eommunieation 
means  ineluding  direet  links  (e.g.,  gestures,  optieal,  sateom,  reg¬ 
ular  phone  eall,  e-mail),  soeial  media  (e.g.,  Faeebook,  Twitter, 
blogs),  mailing  lists,  and  newsgroups.  Sueh  a  network  may  be 
modeled  as  a  hypergraph  with  the  weight  of  a  hyperedge  re¬ 
fleeting  the  eost,  eredibility,  and/or  delay  in  disseminating  infor¬ 
mation  among  all  vertiees  of  this  hyperedge.  In  partieular,  the 
weight  of  a  hyperedge  ean  eapture  the  unique  effeet  on  the  infor¬ 
mation  after  it  passes  through  a  group  of  people.  For  instanee, 
a  result  ean  be  diseussed  by  overlapping  blog  eollaboration  net¬ 
works  as  it  spreads,  and  often  the  diseussion  yields  a  better  result 
than  if  it  only  spreads  through  individuals.  The  minimum-eost 
information-passing  in  soeial  networks  ean  thus  be  modeled  as 
a  shortest  hyperpath  problem  with  suitably  defined  weights. 

Another  potential  applieation  is  that  of  finding  the  most  im¬ 
portant  aetor  and  the  latent  relationships  in  a  soeial  network. 
Under  a  graph  model  of  soeial  networks,  the  relative  importanee 
of  a  vertex  ean  be  measured  by  its  betweenness  and  eloseness 
eentrality  indiees.  The  former  is  defined  based  on  the  number 
of  shortest  paths  that  pass  through  this  vertex,  and  the  latter,  the 
total  weight  of  the  shortest  paths  from  this  vertex  to  all  the  other 
vertiees  [13].  In  a  soeial  network  exhibiting  hyper-relationships, 
betweenness  and  eloseness  eentrality,  based  on  the  shortest  hy¬ 
perpaths,  would  be  better  indieators  of  the  relative  importanee 
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of  each  actor.  In  Section  VII,  we  apply  the  proposed  shortest  hy¬ 
perpath  algorithms  to  the  Enron  e-mail  data  set.  We  propose  a 
weight  function  that  leads  to  the  successful  identification  of  the 
CEO  of  Enron  as  the  most  important  actor  under  the  closeness 
centrality  metric.  The  distance  of  each  person  in  the  data  set  to 
the  CEO  along  the  resulting  shortest  hyperpaths  closely  reflects 
the  hierarchical  structure  of  the  company. 

II.  Dynamic  Shortest  Hyperpath  Problem 

In  this  seciton,  we  introduce  some  basic  concepts  of  hyper¬ 
graph  and  define  the  static  and  the  dynamic  shortest  hyperpath 
problems.  Some  basic  properties  of  the  shortest  hyperpaths  are 
established  and  will  be  used  in  developing  the  dynamic  algo¬ 
rithms  in  subsequent  sections. 

A.  Hypergraph  and  Hyperpath 

Let  V  be  a  finite  set  and  E  a  family  of  subsets  of  V .  If  for  all 
elements  Ci  e  E,  the  following  conditions  are  satisfied: 

^  0,  ei  =  V 

then  the  couple  H  =  (V^E)  is  called  a  (undirected)  hypergraph. 
Each  element  v  e  V  is  called  a  vertex,  and  each  element  e  e  E 
a  hyperedge. 

A  weighted  undirected  hypergraph  is  a  triple  H  =  {V,  E,w) 
with  w  :  E  ^  C  {0}}  being  a  nonnegative  weight  function 
defined  for  each  hyperedge  in  E. 

In  a  hypergraph,  a  hyperpath  is  defined  as  follows. 

Definition  1:  A  hyperpath  between  two  vertices  u  and  v  is 
a  sequence  of  hyperedges  {eo,  ei, . . . ,  e^}  such  that  u  G  eo, 
V  G  Cm,  and  ei  fl  /  0  for  z  =  0, . . . ,  m  —  1.  A  hyperpath 
is  simple  if  nonadjacent  hyperedges  in  the  path  are  nonoverlap¬ 
ping,  i.e.,  Ci  n  Cj  =  0,  Vj  ^  {z,z  ±  1}. 

Let  Le  =  {eo, . . . ,  e^}  be  a  hyperpath  in  a  weighted  hyper¬ 
graph  H.  We  define  the  weight  of  Lg  as 

m 

w{Le)  =  V  w(ei). 
i=o 

B.  Shortest  Hyperpath  and  Relationship  Tree 

Given  two  vertices  u  and  v,  a  natural  question  is  to  And  the 
shortest  hyperpath  (in  terms  of  the  path  weight)  from  u  to  v. 
Since  the  weight  function  is  nonnegative,  it  suffices  to  consider 
only  simple  hyperpaths.  If  the  shortest  hyperpath  is  not  simple, 
we  can  always  generate  a  simple  hyperpath  without  increasing 
the  weight  by  deleting  all  the  hyperedges  between  two  overlap¬ 
ping  nonadjacent  hyperedges. 

The  dynamic  shortest  hyperpath  problem  can  be  similarly 
defined  for  a  sequence  C  =  {(^i,  ^2,  •  •  • ,  of  hyperedge 
changes.  Hyperedge  changes  are  of  the  same  four  types  as 
edge  changes  in  a  graph:  weight  increase,  weight  decrease, 
hyperedge  insertion,  and  hyperedge  deletion.  Similarly,  weight 
increase  and  hyperedge  deletion  will  be  treated  together,  so  will 
weight  decrease  and  hyperedge  insertion. 

In  this  paper,  we  consider  the  single-source  shortest  hy¬ 
perpath  problem:  Find  the  shortest  hyperpaths  from  a  given 
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Fig.  2.  Hyperpaths  and  the  associated  relationship  trees. 


source  s  to  all  other  vertices.  We  focus  on  undirected  hyper¬ 
graphs  first  and  defer  the  treatment  of  directed  hypergraphs  to 
Section  VI. 

We  establish  a  basic  property  of  shortest  hyperpaths.  It  states 
that  any  segment  of  a  shortest  hyperpath  is  also  a  shortest  hy¬ 
perpath  between  any  two  vertices  with  each  in  one  of  the  two 
end  hyperedges  of  this  segment. 

Lemma  1:  Let  L  =  {ei,  62, . . . ,  e/}  be  a  shortest  hyperpath 
from  5  G  ei  to  ^  G  e/.  Then,  for  any  vertex  z;  G  D  e^+i,  the 
hyperpath  =  {ei,  62, . . . ,  e^}  is  a  shortest  hyperpath  from  s 
to  V.  Furthermore,  for  any  two  vertices  u,v  G  H  (if  there 
exist  at  least  two  vertices  in  H  e^+i),  the  shortest  distances 
from  5  to  both  vertices  are  equal. 

Proof:  We  will  prove  by  contradiction.  Assume  that 
Ly  =  {ei,  62, . . . ,  e^}  is  not  a  shortest  hyperpath  for  v.  Then, 
there  exists  a  different  hyperpath  =  {e{,  63, . . . ,  e'^} 

with  w{Vf)  <  w{Ly).  Then,  consider  the  hyperpath  V  = 
{ei,e2,...,e'j.,ei+i,ei+2,...,eJ,  we  have  w(L')  <  w{L), 
which  contradicts  the  fact  that  L  is  a  shortest  hyperpath  to 

This  completes  the  proof  for  the  first  part  of  the  lemma. 
Furthermore,  for  any  two  nodes  u^v  G  D  e^+i,  since  Ly  is 
the  shortest  hyperpath  for  both  vertices,  the  shortest  distances 
from  s  to  both  vertices  equal  to  w{Ly).  ■ 

Next,  we  introduce  the  concept  of  relationship  tree  that  is 
needed  in  the  proposed  dynamic  shortest  hyperpath  algorithm 
HE-DSP.  Since  two  adjacent  hyperedges  in  a  hyperpath  may 
overlap  at  more  than  one  vertex,  the  shortest  hyperpaths  from 
s  to  all  other  vertices  do  not  generally  form  a  tree  in  the  orig¬ 
inal  graph  sense.  For  the  development  of  the  dynamic  shortest 
hyperpath  algorithms,  we  introduce  the  concept  of  relationship 
tree  to  indicate  the  parent-child  relationship  along  shortest  hy¬ 
perpaths.  The  concept  can  be  easily  explained  in  the  example 
given  in  Fig.  2.  Let  {ei,  62}  be  a  shortest  hyperpath  from  5  to 
z;4.  By  Lemma  1,  {ei}  is  a  shortest  hyperpath  for  both  vi  and 
V2.  As  illustrated  in  Fig.  2,  there  are  four  possible  relationship 
trees  to  indicate  the  parent-child  relationship  in  these  shortest 
hyperpaths.  We  will  show  in  Section  III  that  the  choice  of  the 
relationship  tree  does  not  affect  the  correctness  or  performance 
of  the  proposed  algorithm  HE-DSP. 

C.  Notations 

In  the  following,  D  [z;]  denotes  the  distance  of  a  vertex  v  to  the 
source  s  on  the  shortest  hyperpath,  P[z;]  the  parent  of  v  in  the 
chosen  relationship  tree  associated  with  the  shortest  hyperpaths. 
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and  E[v]  the  hyperedge  eontaining  v  and  P[v]  on  the  shortest 
hyperpath  (i.e.,  the  hyperedge  that  leads  to  v  from  P[v]  on  the 
shortest  hyperpath). 

A  vertex  v  is  ealled  an  affeeted  vertex  if  D[v],  P[v],  and/or 
E[v]  ehange  in  the  new  set  of  shortest  hyperpaths.  A  hyperedge 
is  ealled  an  affeeted  edge  if  it  eontains  an  affeeted  vertex.  When 
it  is  neeessary  to  distinguish  the  shortest  distanee  before  and 
after  a  weight  ehange,  d[v]  will  denote  the  shortest  distanee  be¬ 
fore  the  ehange,  d'  [?;]  the  shortest  distanee  after  the  ehange,  and 
D  [tj]  the  aetual  value  stored  in  the  data  strueture  during  the  ex- 
eeution  of  the  algorithm. 

III.  Hyperedge-Based  Dynamic  Shortest  Path  (HE-DSP) 
Algorithm 

In  this  seetion,  we  propose  HE-DSP.  It  is  an  extension  of  the 
dynamie  Dijkstra’s  algorithm  to  hypergraphs.  The  extension  is 
more  eomplex  than  Gallo’s  extension  of  the  statie  Dijkstra’s  al¬ 
gorithm  sinee  the  dynamie  Dijkstra’s  algorithm  for  graphs  relies 
on  the  tree  strueture  of  the  shortest  paths,  a  strueture  no  longer 
there  for  the  shortest  hyperpaths. 

A.  Hyperedge  Weight  Decrease 

For  weight  deerease,  eonsider  that  the  weight  of  a  hyperedge 
e  deereases  to  u^new  It  can  be  shown  that  the  vertex  u  e  e 
with  D  [u]  =  min^;^e  ['^] }  will  not  be  affeeted.  We  then  eheek 
whether  the  other  vertiees  in  e  are  affeeted  by  eheeking  whether 
the  following  inequality  holds: 

D[u]  +  <  D[v].  (1) 

We  then  put  all  the  affeeted  vertiees  into  a  priority  queued  Q. 
The  rest  of  the  proeedure  is  similar  to  the  dynamie  Dijkstra’s 
algorithm,  exeept  that  when  we  update  the  distanee  of  a  vertex, 
we  eheek  all  the  hyperedges  that  eontain  this  vertex.  A  detailed 
implementation  of  the  algorithm  is  given  as  follows. 

HE-DSP:  Weight  Decrease  (e,  iPnew)- 
StepO  (Update  the  hypergraph) 

1  uj(e)  •<  rCnew 

Stepl  (Determine  the  affected  vertices  in  e) 

2  X  ^  avgmm^^^{D[v]} 

3  for  each  v  e  e  sueh  that  D  [x]  -h  Wnew  <  D  [t’]  do 

4  D[v]  ^  D[x]  -h  UJnew;  ^  X]  E[v]  ^  6 

5  Enqueue ( (3,  (t’,  D[t’])) 

6  end 

Step2  (Iteratively  enqueue  and  update  affected 
vertices) 

7  while  NonEmpty  (Q)  do 

8  {z^  D[z])  ^Dequene{Q) 

9  for  each  e  e  E  s.t.  zee 

10  for  each  v  e  e 

11  if  D[v]  >  D[z]  -h  w{e)  then 

1 A  priority  queue  is  an  abstract  data  type  with  the  following  aeeess  protoeol: 
Only  the  highest-priority  element  can  be  aecessed.  Basie  operations  of  a  pri¬ 
ority  queue  include  Enqueue  (add  a  new  item  to  the  queue),  Dequeue  (remove 
the  item  with  the  highest  priority  and  return  this  item).  Update  (change  the  pri¬ 
ority  of  one  item  in  the  queue),  and  Peek  (obtain  the  value  of  the  item  with  the 
highest  priority).  Standard  implementations  of  a  priority  queue  with  different 
time  eomplexities  include  array,  link  list.  Binary  heap,  and  Fibonaeci  heap  [14]. 


12  D[v]^D{z)P 

w{e)]  P[v]  ^  z;  E[v]  ^  e 

13  Enqueue  or  Update((3,  {v.  D[v])) 

14  end; end; end; end 

The  following  theorem  states  the  eorreetness  of  the  algorithm. 

Theorem  1:  If  before  the  weight  deerease,  D[v]  =  d[v\,  E[v\, 
and  P[v\  are  eorreet  for  all  v  e  V,  then  after  the  weight  de¬ 
erease,  D[v]  =  d'[v],  and  E[v]  and  P[v]  are  eorreetly  updated 
by  HE-DSP. 

Proof:  See  Appendix  A.  ■ 

B.  Hyperedge  Weight  Increase 

To  handle  hyperedge  weight  inerease,  speeial  eare  must  be 
taken  to  find  all  the  affeeted  vertiees.  This  proeess  in  the  graph 
ease  relies  on  the  tree  strueture  of  the  shortest  paths  that  no 
longer  exists  in  shortest  hyperpaths.  Our  solution  is  to  use  a 
relationship  tree  as  introdueed  in  Seetion  II-B. 

Consider  that  the  weight  of  a  hyperedge  e  inereases  to  u^new  • 
If  e  is  not  an  edge  in  any  of  the  eurrent  shortest  hyperpaths, 
then  none  of  the  vertiees  will  be  affeeted;  all  shortest  hyper¬ 
paths  remain  unehanged.  Otherwise,  the  deseendants,  and  only 
the  deseendants  of  this  edge  in  the  eurrent  relationship  tree, 
may  be  affeeted.  For  these  vertiees,  some  of  them  will  have  in- 
ereased  distanees,  some  of  them  will  go  through  an  alternative 
path  with  the  same  distanee  (but  ehanged  parent  or  hyperedge), 
while  the  rest  will  not  be  affeeted.  In  order  to  elassify  the  ver¬ 
tiees  into  these  three  eategories,  we  propose  the  following  eol- 
oring  proeess,  whieh  is  a  modifieation  of  the  eoloring  idea  in 
Frigioni’s  algorithm  for  graphs  [4]  to  take  eare  of  the  nonunique 
ehoiee  of  parent  in  hypergraphs. 

1)  V  is  eolored  white  if  d'[v]  =  d[v\  while  keeping  the 
eurrent  P[v\  and  E[v]. 

2)  V  is  eolored  pink  if  f[v\  =  d[v\,  but  only  possible 
through  a  new  P[v\  or  E[v\  or  both. 

3)  V  is  eolored  red  if  d'[v]  <  d[v]. 

It  ean  be  shown  that  if  a  vertex  x  is  white  or  pink,  all  its 
deseendants  in  the  relationship  tree  are  white;  if  x  is  red,  all 
its  deseendants  are  either  red  or  pink.  Therefore,  we  have  the 
following  eoloring  proeedure.  We  first  determine  the  eolor  of 
eaeh  vertex  in  e.  Speeifieally,  the  vertex  u  e  e  with  D[u\  = 
mmy^e{D[v]}  will  not  be  affeeted  and  will  be  eolored  white. 
Any  other  vertex  (say  v)  in  e  will  be  either  pink  or  red,  whieh 
we  ean  determine  by  eheeking  whether  there  is  an  alternative 
shortest  hyperpath  with  the  same  distanee  for  v  (note  that  v 
eannot  be  white  due  to  the  weight  ehange  of  hyperedge  e  that  is 
on  its  eurrent  shortest  hyperpath);  if  sueh  a  path  exists,  then  we 
eolor  V  pink,  otherwise  we  eolor  it  red  and  put  all  its  ehildren  in  a 
priority  queue  M.  The  proeedure  then  iterates  over  eaeh  vertex 
in  M  aeeording  to  an  inereasing  order  of  the  vertex  distanees. 

After  the  eoloring  proeess,  we  only  need  to  deal  with  the  red 
vertiees.  For  eaeh  red  vertex  2;,  we  initialize  its  distanee  with  the 
distanee  of  the  shortest  path  through  one  of  its  nonred  neighbors 
and  put  2;  in  another  priority  queue  Q  (if  no  nonred  neighbor 
exists,  we  initialize  it  with  cx)).  We  then  iterate  by  extraeting  the 
vertex  at  the  top  of  Q  and  updating  its  neighbors  and  Q  until  Q 
is  empty. 
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A  detailed  implementation  of  the  algorithm  is  given  as 
follows. 

HE-DSP:  Weight  Increase  (e,  iCnew)- 

StepO  (Update  the  hypergraph) 

1  <  'W^new 

Stepl  (Determine  the  affected  vertices  in  e) 

2  for  each  v  e  e  s.t.  E[v\  —  e  do 

3  Enqueue (M,  (tj, 

Step2  (Coloring  process) 

4  while  NonEmpty(M) 

5  (^,  [2;])  ^  Dequeue (M) 

6  if  3  nonred  q  e  V  s.t.  3e  G  E  with  q,z  e  e 

and  D[q]  w{e)  =  D[z] 

1  then  z  is  pink;  P[z]  =  q;  E[z]  =  e; 

8  else  z  is  red;  Enqueue(M,  all  ^’s  ehildren) 

9  end; end 

Step3.a  (Initialize  the  distance  vector  for  red  vertices) 

10  for  each  red  vertex  ^  do 

11  if  ^  has  no  nonred  neighbor 

12  then  D[z]  < — hcx);  P[z]  ^  Null 

13  else 

14  let  u  be  the  best  nonred  neighbor  of  z 

15  £[2]  ^argmingg£;^g3„_Uw(e)}; 

16  D[z\*- D[u]+w{E[z\)]P[z\^u-, 

17  Enqueue((5,  (^,  D[^])) 

18  end; end;  end 

Step3.b:  Step2  of  HE-DSP:  Weight  Decrease 

The  theorem  below  states  the  eorreetness  of  the  algorithm. 

Theorem  2:  If  before  the  weight  inerease,  D[v]  =  d[v\,  E[v\ 
and  P[u]  are  eorreet  for  all  u  G  V,  then  after  the  weight  inerease, 
D[v]  =  d^[v],  and  also  E[v]  and  P[v]  are  eorreetly  updated  by 
HE-DSP  regardless  of  the  ehoiee  of  the  relationship  tree. 

Proof:  See  Appendix  B.  ■ 

IV.  Dimension-Reduction-Based  Dynamic  Shortest 
Path  (DR-DSP)  Algorithm 

In  this  seetion,  we  propose  DR-DSP.  DR-DSP  reduees  the 
single  hyperedge  ehange  problem  in  a  hypergraph  to  a  bateh 
problem  in  an  underlying  graph.  When  the  dynamie  problem 
degenerates  to  the  statie  problem,  DR-DSP  leads  to  an  alterna¬ 
tive  algorithm  for  solving  the  statie  shortest  hyperpath  problem. 

A.  Static  Case:  DR-SP 

We  first  eonsider  the  statie  version  of  the  algorithm  (referred 
to  as  DR-SP),  whieh  eaptures  the  basie  idea  of  dimension  re- 
duetion.  We  introduee  an  underlying  graph  that  eaptures  the  es¬ 
sential  (rather  than  eomplete)  information  for  shortest  path  eal- 
eulation  with  only  n  vertiees  (see  Seetion  VIII  for  a  eomparison 
of  our  approaeh  and  other  graph  transformation  sehemes). 

The  proposed  DR-SP  algorithm  is  based  on  the  following  the¬ 
orem  in  whieh  we  show  that  for  a  general  hypergraph  H,  the 
weight  C(;(L* )  of  the  shortest  path  L*  of  H  is  equal  to  the  weight 
of  the  shortest  path  Lq  of  an  indueed  graph  G  derived  from  H. 
Speeifieally,  eorresponding  to  every  hyperedge  e  in  H,  G  eon- 
tains  a  elique  defined  on  the  vertiees  of  e. 
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Theorem  3:  Let  H  =  ( V,  E,  u;)  be  a  hypergraph,  and  G  = 
( V,  E)  the  indueed  graph  of  H  where  an  edge  e  G  E  if  and  only 
if  3e  G  E  sueh  that  e  C  e.  For  eaeh  edge  e  in  G,  its  weight 
WG{e)  is  defined  as 

wg{g)  =  min  _  w{e).  (2) 

{e^E:eDe} 

Let  L*  and  Lq  be  the  shortest  paths  from  ueV  to  veV  in  H 
and  G,  respeetively.  Then,  we  have  that 

w{L^)  =  WG{Lh)^ 

Proof:  First,  for  eaeh  shortest  path  Lq  in  G,  we  ean  obtain 
a  eorresponding  hyperpath  L  in  H  with  the  same  weight  based 
on  (2),  therefore  we  have  that 

WG{Lh)  =  w{L)>w{L^). 

Then,  it  suffiees  to  show  that  there  exists  a  path  Lg 
in  G  sueh  that  wg{Lg)  <  w{L*),  whieh  implies  that 

wg{Lg)  ^  '^g{Lg)  <  'w;(L*). 

Assume  that  L*  =  {eo,  ei , . . . ,  Ck-i}  is  a  shortest  hyperedge 
path  from  vq  to  Vk  in  H  where  vq  G  cq  and  Vk  G  Ck-i.  Let 
Vi  G  Ci-i  n  (z  =  1,  2, . . . ,  ^  —  1)  be  one  of  the  vertiees 
in  the  interseetion  of  hyperedges  e^-i  and  e^.  Constmet  a  path 
Lg  =  {uo,  ui, . . . ,  Ufc}  in  the  graph  G.  For  eaeh  edge  ei  = 
{vi^Vi^i}  (z  =  0, 1, . . . ,  A:  —  1),  sinee  C  e^,  it  follows  from 
(2)  that 

wg{G)  < 

Thus 

fc-i  fc-i 

wg{Lg)  =  ^  WGiei)  <  ^  w{ei)  =  u;(L*) 
i=0  i=0 

i.e.,  wg{Lg)  <  ■ 

It  follows  from  Theorem  3  that  the  shortest  path  in  a  general 
hypergraph  ean  be  obtained  by  applying  Dijkstra’s  algorithm  to 
the  indueed  graph  G  as  defined  in  the  theorem. 

B.  Dynamic  Case:  DR-DSP 

We  now  eonsider  the  dynamie  ease  of  the  problem.  Sinee 
the  underlying  graph  introdueed  in  Seetion  IV-A  does  not  eap- 
ture  the  eomplete  information  of  the  original  hypergraph,  the 
key  ehallenge  here  is  how  to  update  the  underlying  graph  with 
polynomial  time  eomplexity  when  ehanges  oeeur  in  the  original 
hypergraph.  Our  approaeh  is  to  introduee  a  speeial  data  stme- 
ture — a  list  of  priority  queues — to  maintain  the  eomplete  infor¬ 
mation  of  the  original  hypergraph. 

In  the  dynamie  ease,  a  sequenee  G  =  {(^1,  (^2,  •  •  • ,  of 
hyperedge  ehanges  in  the  hypergraph  H  results  in  a  sequenee 
of  edge  ehanges  in  the  indueed  graph  G.  For  eaeh  hyperedge 
ehange  6i,  DR-DSP  first  updates  the  indueed  graph  G  to  loeate 
all  the  ehanged  edges  eaused  by  6i.  In  the  next  step,  DR-DSP 
updates  the  shortest  path  tree  in  the  indueed  graph  G. 

Consider  first  the  graph  update.  A  ehange  to  a  hyperedge  e 
only  affeets  those  edges  in  G  that  are  subsets  of  e,  i.e.,  a  hy¬ 
peredge  ehange  is  loealized  in  the  indueed  graph  G.  Further¬ 
more,  sinee  the  weight  of  an  edge  in  G  is  the  minimum  weight 
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of  all  hyperedges  eontaining  it,  not  all  edges  in  G  that  are  sub¬ 
sets  of  e  will  ehange  weight.  Based  on  these  observations,  we 
propose  a  speeial  data  stmeture  and  proeedure  for  updating  the 
indueed  graph  G  without  regenerating  the  graph  from  serateh. 
Speeifieally,  at  the  initialization  stage  of  the  algorithm,  a  priority 
queue  Muv  for  eaeh  pair  of  vertiees  {u^v)  in  the  hypergraph  is 
established  to  store  the  weights  of  all  hyperedges  that  eontain 
both  u  and  v.  When  a  ehange  oeeurs  to  hyperedge  e,  all  the  pri¬ 
ority  queues  Muv  assoeiated  with  the  pair  of  vertiees  {u^v)  that 
are  eontained  in  e  are  updated  with  the  new  weight  of  e.  Thus, 
the  top  of  these  priority  queues  always  maintains  the  weight  for 
edge  (li,  v)  in  the  indueed  graph  G  for  eaeh  (w,  v).  A  detailed 
implementation  of  the  proposed  proeedure  is  given  as  follows. 

Graph  Update  ew^ew 

1  for  each  u,  G  e 

2  Update^ ,  <C  G^Wuew 

2  '^uv  ^ 

4  end; 

After  the  indueed  graph  G  has  been  updated,  we  now  faee 
a  dynamie  shortest  path  problem  in  a  graph.  However,  sinee  a 
single  hyperedge  ehange  ean  result  in  multiple  edge  ehanges  in 
G,  we  need  to  handle  a  bateh  problem.  While  existing  bateh  al¬ 
gorithms  and  iterative  single-ehange  algorithms  for  graphs  ean 
be  direetly  applied  here,  we  show  that  the  bateh  problem  we 
have  at  hand  has  two  unique  properties  that  ean  be  exploited  to 
improve  the  effieieney  of  the  algorithm. 

Property  1:  The  edge  ehanges  in  G  eaused  by  a  hyperedge 
ehange  are  either  all  weight  deereases  or  all  weight  inereases. 

Property  2:  All  ehanged  edges  in  G  eaused  by  a  hyperedge 
ehange  belong  to  a  elique  in  G. 

In  the  ease  of  weight  deerease,  if  the  weight  of  hyperedge  e 
deereases  to  u^new  ,  by  Theorem  3  and  Property  1,  there 
are  (possibly)  several  edge-weight  deereases  in  the  indueed 
graph  G.  We  know  that  there  is  at  least  one  unaffeeted  node 
X  =  argmin^^g{D[^J]}.  By  Property  2,  these  affeeted  edges 
are  eontained  in  a  elique  derived  from  the  ehanged  hyperedge; 
therefore,  it  is  sufiheient  to  determine  the  distanee  of  every 
node  V  (other  than  x)  in  the  original  ehanged  hyperedge  e  by 
eheeking  whether  D [x]  -h  w^ew  <  D[v].  We  ean  initialize  the 
priority  queue  with  those  nodes  whose  weights  deerease.  After 
that,  the  proeedure  is  similar  to  dynamie  Dijkstra’s  algorithm  in 
the  graph  ease.  Detailed  implementations  of  DR-DSP  weight 
deerease  are  given  as  follows. 

DR-DSP:  Weight  Decrease  (e,  u^new)- 

StepO  (Update  the  hypergraph  and  G) 

1  w(^e^  <  '^new 

2  Graph  Update  (e,  w^ew) 

Stepl  of  HE-DSP:  Weight  Decrease 

Step2  (Iteratively  update  all  affected  vertices) 

3  while  NonEmpty  (Q)  do 

4  {z^  D[z])  ^Dequene{Q) 

5  for  each  v  e  V  s.t.  {z,  v)  e  E 

6  if  D[v]  >  D[z]  w{z^  v)  then 

7  D[v\^  D{z) -\- w{z^v)]P[v\^  z 

8  Enqueue  or  Update(Q,  ^H)) 
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9  end; end;  end 

In  the  ease  of  weight  ineresae,  if  the  weight  of  hyperedge  e 
inereases  to  ujnew,  by  Theorem  3  and  Property  1,  there  are  (pos¬ 
sibly)  several  edge- weight  inereases  in  the  indueed  graph  G. 
Similar  to  the  single-ehange  ease  in  graph,  there  is  at  least  one 
unaffeeted  node  x  =  Then,  another  node 

V  e  eis  affeeted  only  if  E[v]  =  e,  i.e.,  e  is  on  its  shortest  hyper¬ 
path.  We  use  all  sueh  nodes  to  initialize  the  priority  queue  M. 
The  rest  is  similar  to  that  in  the  dynamie  Dijkstra’s  algorithm 
ease. 

Detailed  implementations  of  DR-DSP  weight  inerease  are 
given  as  follows. 

DR-DSP:  Weight  (e,  Wnew)- 

StepO  (Update  the  hypergraph  and  G) 

1  ic(e)  ■<  U^new 

2  Graph  Update  (e,  u;new) 

Stepl  of  HE-DSP:  Weight  Increase 

Step  2  (Coloring  Process) 

3  while  NonEmpty(M) 

4  (^,  D[^])  ^  Dequeue(M) 

5  if  3  nonred  g  G  U  s.t.  D[q]  +  w{q,  z)  —  D[z] 

6  then  z  is  pink 

7  else  z  is  red;  Enqueue(M,  all  ^’s  ehildren) 

8  end; end 

Step3.a  (Initialize  the  distance  vector  for  red  vertices) 

9  for  each  red  vertex  ^  do 

10  if  ^  has  no  nonred  neighbor 

11  then  D[z]  ^  +cx);  P[z]  ^  Null 

12  else 

13  let  u  be  the  best  nonred  neighbor  of  z 

14  D[z]^  D[u] -\- w{u^  z)]  P[z]^  u 

15  Enqueue((5,  (^,  D[^])) 

16  end; end;  end 

Step3.b:  Step2  of  DR-DSP:  Weight  Decrease 


V.  Time  Complexity  Analysis 

Given  a  hypergraph  H  =  (V^E^w)  and  a  ehange  to  hyper¬ 
edge  e,  let  \6\  denote  the  number  of  affeeted  vertiees,  \6^\  = 
EeGis.e  is  affected  kP,  m  the  number  of  edges  in  the  induced 
graph,  and  ||^||  the  number  of  affeeted  edges  in  the  indueed 
graph  plus  \6\. 

Theorem  4:  The  time  eomplexities  of  HE-DSP  and  DR-DSP 
are  0(|^|  log  |^|  -h  |(5^|)  and  0(|^|  log \6\  +  ||^||  -h  |ep  logm), 
respeetively. 

Proof:  The  time  eomplexity  of  DR-SP  mainly  eomes  from 
Steps  1  and  2.  Step  2  is  essentially  applying  Dijkstra’s  algo¬ 
rithm  to  a  graph  with  n  vertiees  and  m  edges  where  m  is  the 
number  of  edges  in  the  indueed  graph  G.  The  running  time  is 
thus  0(n  log n  -h  m).  An  implementation  of  Step  1  is  to  ob¬ 
tain  the  edge  weight  wg{E)  based  on  (2).  Therefore,  the  time 
eomplexity  for  Step  1  is  kP)?  ^(^)-  ^bh  m 

upper-bounded  by  (sinee  for  eaeh  e  G  E,  there  are  at  most 
|e|(|e|  —  l)/2  edges  in  G),  we  arrive  at  the  total  time  eomplexity 
ofDRSP. 


This  article  has  been  accepted  for  inclusion  in  a  future  issue  of  this  journal.  Content  is  final  as  presented,  with  the  exception  of  pagination. 


GAO  et  al.\  DYNAMIC  SHORTEST  PATH  ALGORITHMS  FOR  HYPERGRAPHS 


7 


For  Gallo’s  Algorithm,  similar  to  Dijkstra’s  algorithm,  the 
time  eomplexity  is  mainly  in  updating  the  neighbors  of  the  non- 
fixed  vertex  2:  with  the  minimal  distanee  D[z\.  For  eaeh  2;,  the 
algorithm  seans  all  the  hyperedges  eontaining  2;.  For  eaeh  pair 
of  vertiees  (u,  v)  £  e,  e  is  seanned  twiee.  Therefore,  the  total 
number  of  sueh  operations  is  =  XlecE;  |ep.  Also,  extraeting  2: 
from  the  priority  queue  implemented  by  a  Fibonaeei  heap  takes 
O(logn)  time.  The  total  time  eomplexity  of  Gallo’s  algorithm 
thus  follows. 

For  a  simplieial  eomplex,  =  0{d?2'^m),  the  eomplexity 
of  Gallo’s  algorithm  thus  follows.  For  DR-SP,  exploiting  the 
property  that  the  edge  set  is  elosed  under  the  subset  operation  in 
a  simplieial  eomplex,  we  ean  use  a  top-down  seheme  in  Step  1  of 
DR-SP  to  ealeulate  the  weight  wg{s)  induetively  with  respeet 
to  the  dimension  of  a  faeet  as  follows: 

'^g{s)  —  min{r(;(5),  s  D  5,  dim[5']  =  z  +  1}} 

where  wq  {s')  =  w {s')  for  the  faeet  s' .  The  time  eomplexity  for 
Step  1  ean  then  be  improved.  Beeause  eaeh  z-dimensional  faee 
is  assoeiated  with  d  —  i  eomparisons,  the  running  time  of  Step  1 
for  eaeh  ti-dimensional  faeet  is  given  by 


Therefore,  the  time  eomplexity  for  Step  1  is  0{d2'^m) .  The  total 
time  eomplexity  thus  follows.  ■ 

From  Theorem  4,  we  see  that  if  6  is  small  and  |e|  is  large, 
HE-DSP  performs  better,  sinee  in  DR-DSP,  the  indueed  graph 
has  to  be  updated  regardless  of  whether  there  are  affeeted  ver¬ 
tiees.  Thus,  in  a  sequenee  of  hyperedge  ehanges,  if  only  a  small 
fraetion  of  them  aetually  have  affeeted  nodes,  then  HE-DSP  will 
outperform  DR-DSP.  On  the  other  hind,  if  8  is  large,  DR-DSP 
will  outperform  HE-DSP  sinee  usually  >>  ||^||. 

VL  Directed  Hypergrpah 

In  this  seetion,  we  extend  both  algorithms  to  direeted 
hypergraphs. 


B.  Extensions 

First,  we  extend  HE-DSP  to  direeted  hypergraphs.  Consider 
first  the  ease  with  weight  deerease.  In  Step  1,  we  still  identify 
the  affeeted  vertiees  in  the  ehanged  hyperedge  e.  However,  sinee 
the  ehanged  hyperedge  has  only  one  souree  vertex,  if  any  other 
vertex  is  affeeted,  its  new  shortest  path  has  to  eome  from  the 
souree  vertex.  Speeifieally,  if  a  vertex  u  satisfies  D[Se]  + 
z^new  <  D[u],  then  it  has  a  shortest  path  from  Se  through  the 
ehanged  hyperedge  e  and  beeomes  an  affeeted  vertex.  We  then 
put  it  into  the  priority  queue  Q.  In  Step  2,  we  still  extraet  the 
vertex  2;  with  minimal  eurrent  distanee  in  Q  and  use  it  to  update 
neighbors.  Sinee  the  model  is  direeted,  instead  of  looking  at 
all  hyperedges  that  eontain  2;,  we  only  eonsider  those  direeted 
hyperedges  soureed  at  2;.  The  rest  of  the  algorithm  remains  the 
same. 

For  weight  inerease,  the  extension  of  HE-DSP  is  similar  to 
the  above  proeess.  In  Step  1,  we  put  all  the  reeeiving  vertiees  in 
M  to  initialize  the  eoloring.  Then,  in  both  eoloring  and  updating 
proeess,  when  eheeking  or  updating  the  state  of  neighbors,  we 
simply  ehange  the  set  of  neighbors  through  one  hyperedge  e  to 
the  reeeiving  set  Te . 

Next,  we  eonsider  an  extension  of  the  DR-DSP  algorithm. 
Sinee  the  original  hypergraph  H  is  direeted,  the  indueed  graph 
G  beeomes  a  direeted  graph.  The  weight  of  one  edge  u  ^  v  is 
the  minimum  weight  among  all  hyperedges  that  provide  a  path 
from  uto  V.  More  speeifieally,  the  weight  Wu,v  is  given  by 


'^u,v 


min  w{e). 
6  C  E  '.U^  S  g  j'UCT'e 


(3) 


A.  Definitions 

In  a  direeted  hypergraph  H  =  {V,E),  eaeh  hyperedge  e  e  E 
eonsists  of  a  souree  vertex  Se  and  a  set  of  terminal  vertiees  Te 
(see  Fig.  3  for  an  illustration).  A  direeted  hyperpath  L  from  5  to  t 
is  a  sequenee  of  direeted  hyperedges  L  =  {ei,  62, . . . ,  e^}  sueh 
that  Se^  =  s,  t  E  Tej^  and  Se^  C  Te-_-^  for  z  =  2, . . . ,  A:.  See 
Fig.  3  for  an  example  of  a  hyperpath  with  two  hyperedges:  L  = 
{61,62}  where  S'ei  =  l.T^^  =  (2, 3, 4, 5),  =  3,Te2  = 

{4,6,7}. 

Based  on  the  definition,  one  ean  easily  see  that  an  undireeted 
hypergraph  ean  be  seen  as  a  speeial  ease  of  a  direeted  hyper¬ 
graph.  Speeifieally,  eaeh  undireeted  hyperedge  e  with  k  vertiees 
ean  be  treated  as  k  direeted  hyperedges,  eaeh  from  a  souree 
vertex  z;  in  e  to  the  destination  set  e  —  z;.  Additionally,  these 
k  direeted  hyperedges  share  the  same  weight  as  the  undireeted 
hyperedge. 


Theorem  3  ean  be  easily  extended  to  the  direeted  ease.  Henee, 
the  problem  is  redueed  to  a  speeial  bateh  dynamie  shortest  path 
problem  in  direeted  graph,  and  both  DR-DSP:  Weight  Deerease 
and  DR-DSP:  Weight  Inerease  ean  be  extended. 

VII.  Simulation  Examples 

We  present  simulation  results  on  the  running  time  of  the  pro¬ 
posed  dynamie  shortest  hyperpath  algorithms.  We  test  the  pro¬ 
posed  algorithms  on  hypergraphs  generated  from  a  random  geo- 
metrie  model  as  well  as  those  generated  by  the  Enron  e-mail  data 
set.  All  simulation  eode  is  eompiled  and  run  on  the  same  laptop 
equipped  with  a  3.0-GHz  i7-920XM  Mobile  Proeessor. 

A.  Random  Geometric  Hypergraph 

We  first  eonsider  a  random  geometrie  hypergraph  model  in 
whieh  n  nodes  are  uniformly  distributed  in  an  a  x  a  square.  All 
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nodes  within  a  circle  with  radius  r  form  a  hyperedge  (circles 
are  centered  on  an  x  grid).  The  weight  of  each  hyperedge 
is  given  by  the  average  distance  between  all  pairs  of  vertices  of 
this  hyperedge. 

A  sequence  of  changes  is  then  generated,  and  the  proposed 
dynamic  algorithms  are  employed  to  maintain  all  the  shortest 
hyperpaths  from  the  source  s  located  at  a  comer  of  the  a  x  a 
square.  Each  change  can  be  a  hyperedge  insertion  (with  prob¬ 
ability  pi),  a  hyperedge  deletion  (with  probability  pd),  or  a 
weight  change  (with  probability  1  —  —  po)  with  new  weight 

chosen  uniformly  in  w;max]-  In  the  case  of  a  hyperedge 

deletion  or  a  weight  change,  the  hyperedge  to  be  deleted  or  to  be 
assigned  a  new  weight  is  chosen  according  to  the  two  models  de¬ 
tailed  below.  Hyperedge  insertions  are  only  realized  when  there 
are  hyperedges  that  have  been  deleted,  and  a  randomly  chosen 
one  is  inserted  back.  This  ensures  that  all  hyperedges  satisfy  the 
geometric  property  determined  by  r  at  all  time.  It  also  models 
the  practical  scenario  where  a  broken  link  is  repaired. 

In  selecting  a  hyperedge  for  deletion  or  weight  change,  we 
consider  two  different  models:  the  random  change  model  and 
the  targeted  change  model.  In  the  former,  the  hyperedge  is  ran¬ 
domly  and  uniformly  chosen  among  all  hyperedges.  In  the  latter, 
it  is  randomly  and  uniformly  chosen  from  the  current  shortest 
hyperpaths.  This  models  the  scenarios  where  hyperedges  in  the 
shortest  hyperpaths  are  more  prone  to  changes  due  to  attacks, 
frequent  use,  or  higher  priority  in  maintenance  and  upgrade. 

In  Fig.  4,  we  show  the  simulation  results  on  the  mnning 
time  of  the  two  proposed  algorithms  under  a  sequence  of 
lO'^  changes.  We  see  that  HE-DSP  has  lower  complexity  in 
networks  with  random  topological  and  weight  changes  (Fig.  4, 
top),  whereas  DR-DSP  should  be  preferred  in  networks  with 
targeted  changes  (Fig.  4,  bottom).  This  partition  of  the  appli¬ 
cation  space  can  be  explained  from  the  structures  of  these  two 
algorithms.  Under  the  random  change  model,  a  large  fraction 
of  changes  does  not  result  in  changes  in  the  current  shortest 
hyperpaths.  Such  changes  lead  to  little  computation  in  main¬ 
taining  the  shortest  hyperpaths  for  both  algorithms,  but  require 
about  the  same  amount  of  computation  in  the  Graph-Update 
step  of  DR-DSP  for  maintaining  the  induced  graph.  On  the 
other  hand,  under  the  targeted  change  model,  all  hyperedge 
deletions  and  weight  changes  affect  the  shortest  hyperpaths. 
Updating  the  shortest  hyperpaths  can  be  done  more  efficiently 
in  DR-DSP  since  it  works  on  the  induced  graph  with  a  much 
smaller  number  of  edges. 

B.  Energy-Efficient  Routing  in  Multichannel  Multiradio 
Networks 

Consider  a  multichannel  multiradio  ad  hoc  network  with 
n  nodes  and  m  channels.  Each  node  i  is  associated  with  a 
fixed  transmission  range  Ri,  a  transmission  energy  cost 
a  receiving  energy  cost  and  a  list  of  <  m  channels 
{Cfci,  Cfc2, . . . ,  C  {ci,  C2, ....  Cm}  available  to  it. 

This  network  can  be  modeled  by  a  directed  hypergraph  where 
a  directed  hyperedge  e  consists  of  one  source  vertex  Se  and 
a  set  of  destination  vertices  Tg.  Specifically,  each  node  i  is  a 
vertex  Vi  in  the  hypergraph.  For  each  channel  Cj  available  to 
node  Vi,  there  is  a  hyperedge  Cij  from  Vi  to  the  set  Vij  of  vertices 
consisting  of  nodes  that  are  within  the  transmission  range  of  Vi 


Fig.  4.  Average  mnning  time,  (top)  Random  change  model,  (bottom)  Targeted 
change  model  {n  —  1000,  a  =  1000,  r  =  \/l000,  h  =  1,  pi  = 

Wmin  =  10,  rcmax  =  20,  thc  avcragc  is  taken  over  50  random  hypergraphs). 


and  share  the  channel  Cj  (if  Vij  is  empty,  we  do  not  consider  this 
hyperedge).  Then,  the  cost  of  this  hyperedge  is  defined  as 

w{eij)  =  ^  Pj.  (4) 

veVij 

We  consider  a  mobile  network.  In  this  case,  all  four  types  of 
hyperedge  changes  (i.e.,  edge  insertion,  edge  deletion,  weight 
increase,  and  weight  decrease)  can  occur,  and  multiple  changes 
can  occur  simultaneously.  We  thus  have  a  fiilly  dynamic  batch 
problem  (as  stated  earlier,  the  proposed  algorithms  can  handle 
batch  changes  with  minor  changes  in  their  implementations). 

In  the  simulation  results  shown  in  Fig.  5,  nodes  are  uni¬ 
formly  distributed  within  an  a  x  a  square.  The  list  of  channels 
available  to  each  node  is  drawn  from  all  channels  with  proba¬ 
bility  p.  At  each  time,  one  node  is  randomly  chosen  to  move. 
The  movement  (Aa;,  Ay)  is  generated  by  a  random  vector 
(A A,  Ay).  We  have  implemented  five  algorithms:  repeating 
static  algorithm,  repeating  single  change  HE-DSP,  repeating 
single  change  DR-DSP,  batch  HE-DSP,  and  batch  DR-DSP.  As 
we  can  see  from  Fig.  5,  all  dynamic  algorithms  outperform  the 
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Fig.  5.  Average  running  time  with  parameters  n  =  1000,  m  =  20, 

a  =  R  17(0.8,1.2),  a  ~  17(8,12),  /?  ~  17(0.8,1.2), 

AX,  AY  ~  17(— 0.5,  0.5),  and  p  =  0.2.  The  average  is  taken  over  50 
random  hypergraphs. 


static  algorithm.  Also,  the  batch  algorithms  perform  better  than 
repeating  the  single  change  dynamic  algorithms.  In  both  single 
and  batch  algorithms,  HE-DSP  are  better  than  DR-DSP.  This 
is  mainly  because  many  hyperedge  changes  occur  outside  the 
current  shortest  hyperpaths. 

C.  Enron  E-Mail  Data  Set 

In  this  example,  we  consider  the  application  of  the  shortest 
hyperpath  algorithms  in  finding  the  most  important  actor  in  a 
social  network.  We  consider  the  Enron  e-mail  data  set  and  use 
the  same  hypergraph  generation  model  as  in  [15].  Specifically, 
each  person  is  a  vertex  of  the  hypergraph,  and  the  sender  and 
recipients  of  every  e-mail  form  a  hyperedge.  Our  objective  is 
to  identify  the  most  important  person  measured  by  the  close¬ 
ness  centrality  index  (i.e.,  the  total  weight  of  the  shortest  hyper¬ 
paths  from  this  person  to  all  the  other  persons).  The  first  step 
is  to  assign  weight  to  each  hyperedge  that  reflects  “distance.” 
While  there  is  no  universally  accepted  way  of  measuring  dis¬ 
tance  in  a  social  network  observed  through  e-mail  exchanges, 
certain  general  rules  apply.  First,  a  direct  e-mail  exchange  be¬ 
tween  two  persons  indicates  a  stronger  tie  than  an  e-mail  sent  to 
a  large  group.  Thus,  the  weight  of  a  hyperedge  should  be  an  in¬ 
creasing  function  of  the  cardinality  of  this  hyperedge.  Second, 
more  frequent  e-mail  exchange  among  a  given  group  of  people 
shows  stronger  ties.  Thus,  the  weight  of  a  hyperedge  should  be 
decreasing  with  the  number  of  times  that  this  hyperedge  appears 
in  the  e-mail  data  set.  Considering  these  two  general  rules,  we 
adopt  the  following  weight  function: 

w(e)  =  (\/ie[)“‘'  (5) 

where  |e|  is  the  cardinality  of  the  hyperedge  e,  and  a  is  the 
parameter  measuring  how  fast  the  weight  decreases  with  the 
number  I  of  times  that  this  hyperedge  appears  in  the  data  set. 
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Fig.  6.  Average  distance  from  the  CEO  to  others  at  different  positions. 


We  can  then  apply  DR-SP  on  the  resulting  (static)  hypergraph 
to  find  the  shortest  hyperpaths  rooted  at  each  vertex  and  com¬ 
pute  this  vertex’s  closeness  centrality  index.  With  the  weight 
function  given  in  (5)  using  a  =  0.6,  the  identified  most  impor¬ 
tant  actor  is  the  CEO  of  Enron.  The  average  distance  (along  the 
shortest  hyperpath)  from  the  CEO  to  other  persons  at  various  po¬ 
sitions  is  shown  in  Fig.  6.  We  observe  that,  in  general,  the  higher 
the  position,  the  shorter  the  distance.  These  results  demonstrate 
that  the  adopted  hypergraph  model  and  weight  function  capture 
the  essence  of  the  problem. 

To  demonstrates  advantage  of  the  hypergraph  model,  we  run 
the  same  process  based  on  the  following  graph  model.  Each 
person  is  a  vertex  of  the  graph.  There  exists  an  edge  between 
two  vertices  li,  v  only  if  there  is  an  e-mail  that  involves  the  two 
people  represented  hy  u.  v  either  as  sender  or  receiver.  In  com¬ 
parison  to  the  hypergraph  model,  we  substitute  |e|  in  (5)  to  2, 
the  cardinality  of  an  edge,  in  the  weight  function 

w{u,v)  =  \  (6) 

In  this  model,  the  weight  of  an  edge  only  depends  on  the  fre¬ 
quency  of  e-mail  exchanges,  but  not  the  number  of  people  in¬ 
volved  in  those  e-mails.  First,  we  use  Dijkstra’s  algorithm  to 
compute  closeness  centrality  for  each  vertex.  Based  on  the  re¬ 
sults,  the  most  important  actor  is  still  identified  to  be  the  CEO  of 
Enron.  However,  the  average  distances  from  the  CEO  no  longer 
reveal  the  hierarchy  of  this  company  for  the  entire  range  of  a 
(values  from  0. 1  to  0.9  with  a  step  size  of  0. 1  were  tested).  Fig.  7 
shows  the  result  for  g  =  0.6.  The  result  is  similar  for  other 
values  of  a. 

Next,  we  construct  a  dynamic  hypergraph  sequence  based 
on  the  Enron  data  set.  At  the  beginning,  the  hypergraph  con¬ 
tains  only  individual  vertices.  We  then  consider  each  e-mail 
chronologically.  Each  e-mail  either  adds  a  new  hyperedge  or 
decreases  the  weight  of  an  existing  hyperedge  (due  to  the  in¬ 
creased  number  of  appearances  of  this  hyperedge).  The  two  pro¬ 
posed  algorithms  are  employed  to  maintain  the  shortest  hyper¬ 
paths  rooted  at  the  CEO  after  each  change.  The  running  time  is 
given  in  Fig.  8,  which  shows  the  lower  complexity  of  DR-DSP. 


This  article  has  been  accepted  for  inclusion  in  a  future  issue  of  this  journal.  Content  is  final  as  presented,  with  the  exception  of  pagination. 
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Fig.  9.  Given  a  hypergraph  (a),  we  first  create  a  multigraph  (b)  by  replacing 
each  hyperedge  with  a  clique  whose  edges  have  the  same  weight  as  the  hyper¬ 
edge.  Then,  each  edge  in  this  multigraph  is  replaced  with  one  vertex  and  two 
edges  to  create  a  traditional  graph  (c).  The  weight  of  a  new  edge  in  the  graph  is 
half  the  weight  of  the  multigraph  edge  from  which  the  new  edge  is  derived.  As 
can  be  seen,  the  shortest  path  in  the  induced  graph  represents  the  shortest  path 
in  the  original  hypergraph.  Maintaining  this  underlying  graph  after  one  hyper¬ 
edge  change  takes  )  time.  However,  the  number  of  vertices  in  the  induced 
graph  is  more  than  the  sum  of  the  number  of  vertices  and  the  number  of  hyper¬ 
edges  in  the  original  hypergraph. 


Fig.  7.  Average  distance  from  the  CEO  to  others  at  different  positions  obtained 
in  the  graph  model  (a  =  0.6). 


number  of  hyperedge  changes 


X  10^^ 


Fig.  8.  Average  running  time  for  the  Enron  data  set  (ck  =  0.6,  the  average  is 
taken  over  50  Monte  Carlo  runs). 
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Fig.  10.  Given  a  hypergraph  H  (a),  we  create  its  dual-transformation  bipartite- 
graph  G  (b):  Each  vertex  on  the  left  side  of  G  represents  a  vertex  in  the  H , 
and  each  vertex  on  the  right  side  represents  a  hyperedge  in  H .  There  is  an  edge 
between  two  vertices  if  the  corresponding  hyperedge  contains  the  corresponding 
vertex  in  H .  The  weight  of  the  edge  in  G  is  half  the  weight  of  the  hyperedge 
represented  by  the  right  vertex  of  this  edge.  It  is  easy  to  see  that  the  shortest 
path  problem  in  H  is  equivalent  to  the  shortest  path  in  G.  Maintaining  this 
underlying  graph  after  one  hyperedge  change  takes  0{n)  time.  However,  the 
number  of  vertices  in  the  induced  graph  equals  the  sum  of  the  number  of  vertices 
and  the  number  of  hyperedges  in  the  original  hypergraph. 


The  reason  is  that  a  large  fraetion  of  hyperedge  ehanges  result 
in  ehanges  in  the  shortest  hyperpaths. 

VIIL  Discussion 

The  shortest  path  problem  in  a  hypergraph  ean  be  transformed 
to  the  shortest  path  problem  in  an  indueed  graph,  and  there 
are  many  ways  of  transforming  the  problem.  As  an  algorithmie 
study  of  the  problem  where  the  objeetive  is  the  effieieney  of 
the  solution,  we  aim  to  address  the  following  two  questions  in 
this  paper:  1)  When  is  it  more  efiheient  to  transform  the  problem 
into  a  graph  problem  rather  than  working  direetly  on  the  hyper¬ 
graph?  2)  When  it  is  more  efiheient  to  transfiorm  the  problem, 
whieh  graph  transfiormation  has  better  time  eomplexity? 

The  first  question  was  addressed  by  a  detailed  eomplexity 
eomparison,  both  analytieally  (Seetion  V)  and  empirieally 
(Seetion  VII),  ofi  the  two  proposed  algorithms,  one  working 
direetly  on  the  hypergraph,  the  other  based  on  a  graph  transfior¬ 
mation.  This  eomplexity  eomparison  leads  to  a  partition  ofi  the 


applieation  domain  fior  these  two  approaehes  based  on  the  type 
ofi  the  hypergraph  and  the  nature  ofi  the  dynamies  ofi  the  weight 
ehanges. 

To  taekle  the  seeond  question,  we  first  notiee  that  to  eap- 
ture  all  the  infiormation  ofi  an  edge-weighted  hypergraph  with 
n  vertiees,  the  indueed  graph  ean  be  exponentially  large  with 
n  sinee  the  number  ofi  hyperedges  (henee  hyperedge  weights) 
ean  be  exponential  with  n.  Figs.  9  and  10  demonstrate  two  ex¬ 
amples  ofi  a  graph  transfiormation  that  ean  result  in  an  exponen¬ 
tially  large  graph.  Sinee  time  eomplexity  ofi  dynamie  shortest 
path  algorithms  ean  be  the  same  as  a  statie  one  in  the  worst 
ease  [2],  [4],  [5],  applying  them  to  these  types  ofi  underlying 
graph  that  represents  eomplete  infiormation  ofi  the  given  hyper¬ 
graph  results  in  exponential  time  eomplexity  with  n 

When  adapting  a  graph  reduetion-based  approaeh,  we  do 
need  to  maintain  the  eomplete  infiormation  ofi  the  hypergraph  in 
order  to  update  the  hypergraph  fior  fiuture  hyperedge  ehanges. 
However,  we  also  want  to  reduee  the  size  ofi  the  underlying 
graph  so  that  applying  a  dynamie  shortest  path  algorithm  to  it 
results  in  polynomial  time  eomplexity  with  n  in  the  worst  ease. 


This  article  has  been  accepted  for  inclusion  in  a  future  issue  of  this  journal.  Content  is  final  as  presented,  with  the  exception  of  pagination. 
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TABLE  I 

Comparison  Between  Three  Graph-Transformation-Based 
Approaches 


Multi¬ 
graph 
approach 
in  Fig.  9 

Bipartite 
graph 
approach 
in  Fig.  10 

DR-DSP 

Number  of  Vertices  in  G 

0(2")  1 

0(2") 

n 

Maintaining  G 

O(n^) 

0(n) 

0{n^) 

Applying  DSP  to  G 

0(2") 

0(2") 

O(n^) 

Total  Time  Complexity 

0(2") 

0(2"") 

O(n^) 

Comparison  of  the  three  graph-reduction  based  approaches  while 
there  are  0(2’^)  hyperedges  in  the  hypergraph  and  one  hyperedge 
change  causes  all  vertices  affected.  The  first  row  is  the  number  of 
vertices  in  the  resulting  graph  G  after  the  transformation.  The  second 
row  is  the  time  complexities  for  maintaining  the  underlying  graphs. 
The  third  row  is  the  time  complexities  of  applying  dynamic  shortest 
path  algorithms  to  the  underlying  graph  G.  The  last  row  is  the  total 
time  of  the  corresponding  graph-transformation  based  approach. 

We  overcome  this  dilemma  in  the  proposed  DR-DSP  algo¬ 
rithm  by  constructing  a  two-level  data  structure.  One  level  is  a 
list  of  priority  queues  for  maintaining  the  complete  information 
of  the  hypergraph.  The  other  level  is  an  underlying  graph  with 
n  vertices  that  only  contains  the  essential  information  needed 
for  the  shortest  path  calculation.  Though  the  size  of  the  data 
structure  at  the  first  level  can  be  exponential  with  n,  main¬ 
taining  it  and  the  underlying  graph  after  each  hyperedge  change 
only  requires  polynomial  time  complexity  with  n  (specifically, 
O(n^)  in  the  worst  case).  Since  the  underlying  graph  has  only 
n  vertices,  updating  the  shortest  paths  after  each  change  also 
has  a  time  complexity  polynomial  with  n.  Table  I  gives  a  de¬ 
tailed  complexity  comparison  of  the  proposed  DR-DSP  algo¬ 
rithm  with  the  two  alternative  graph-transformation-based  ap¬ 
proaches  given  in  Figs.  9  and  10.  It  shows  that  DR-DSP  reduces 
the  time  complexity  from  exponential  to  polynomial  by  striking 
a  balance  between  the  complexity  for  maintaining  the  resulting 
graph  G  after  the  transformation  and  the  complexity  for  main¬ 
taining  the  shortest  paths  in  G. 

IX.  Conclusion 

We  have  presented  the  first  study  of  the  fully  dynamic  shortest 
path  problem  in  a  general  hypergraph.  We  have  developed  two 
dynamic  algorithms  for  finding  and  maintaining  the  shortest  hy¬ 
perpaths.  These  two  algorithms  complement  each  other  with 
each  one  preferred  in  different  types  of  hypergraphs  and  net¬ 
work  dynamics,  as  illustrated  in  the  time  complexity  analysis 
and  simulation  experiments.  The  special  two-level  data  struc¬ 
ture  of  our  second  algorithm,  DR-DSP,  may  also  find  other  ap¬ 
plications  in  dynamic  hypergraph  studies.  We  have  discussed 
and  studied  via  experiments  over  a  real  data  set  the  potential 
applications  of  the  dynamic  shortest  hyperpath  problem  in  so¬ 
cial  and  communication  networks. 

Appendix  A 
Proof  of  Theorem  1 

The  proof  is  based  on  the  following  three  lemmas. 
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Lemma  2:  Let  a;  =  argmin^^g{fi[t;]},  then  =  d![x\  and 
d'[x\  =  mmy^e{d'[v]}. 

Proof:  Proof  by  contradiction.  Assume  that  G?' [a;]  <  d[x\, 
then  X  has  to  use  e  on  its  new  shortest  hyperpath.  Since  we 
consider  only  simple  hyperpaths  and  x  e  e,wQ  have  E[x]  = 
e.  Therefore,  its  parent  y  =  P[v]  cannot  use  e  on  its  shortest 
hyperpath,  which  implies  that  the  shortest  distance  to  y  does 
not  change:  d[y]  =  f[y].  Given  that  y  is  the  parent  of  x  on  its 
new  shortest  hyperpath,  we  have  d[y]  =  d'[y]  <  d'[x]  <  d[x], 
which  contradicts  the  definition  of  x. 

For  the  second  statement,  assume  that  there  exists  ^  G  e  such 
that  d!  [^]  <  f  [x] .  Based  on  the  definition  of  x  and  the  hypothet¬ 
ical  assumption,  d[z\  >  d[x]  =  d'[x]  >  d'[z\.  It  thus  follows 
that  ^’s  shortest  hyperpath  changes  and  E[z\  —  e  in  the  new 
shortest  hyperpath.  Following  the  same  line  of  arguments  by 
considering  the  parent  of  we  arrive  at  the  same  contradiction 
in  terms  of  the  definition  of  x.  ■ 

Lemma  3:  For  any  vertex  v,v  is  enqueued  into  Q  if  and  only 
if  d'[v]  <  d[v]. 

Proof:  Consider  first  that  v  is  enqueued  into  Q.  From  the 
algorithm,  this  can  only  happen  if  there  exists  a  neighbor  ^  and 
a  hyperedge  e  3  v^z  such  that  D[z]  +  w{e)  <  D[v].  We  thus 
have  d[v]  >  D[v]  >  D[z]  -h  w{e)  >  D'[v]  (note  that  at  any 
time,  d[v]  >  D[v]  >  D'[v],  which  can  be  easily  seen  from  the 
procedure  of  the  algorithm). 

We  now  prove  the  converse.  Assume  that  f[v]  <  d[v].  Let 
p  =  {ei,  62, . . . ,  ej,  e,  e^+i, . . . ,  ej  be  t;’s  new  shortest  hyper¬ 
path.  There  exists  G  6  n  6^+1  such  that  <  d[ui^i]. 

In  Step  1  of  the  algorithm,  Ui^i  is  enqueued.  Similarly,  there  ex¬ 
ists  Ui^2  ^  ej+i  n  6^+2  with  d'[ui^2]  <  d[ui^2]  -  Then,  Ui^2 
will  be  enqueued  in  Step  2  of  the  algorithm  when  Ui^i  is  de¬ 
queued  if  it  has  not  been  enqueued  before  that.  Repeating  this 
line  of  argument,  we  conclude  that  there  exits  ui  G  e/_i  D  e/ 
with  d'[ui]  <  d[ui]  and  ui  is  enqueued  into  Q.  Then,  v  will  be 
enqueued  when  ui  is  dequeued  if  it  is  not  enqueued  already.  ■ 

Lemma  4:  For  each  v  dequeued  from  Q,  D[v]  =  f[v]. 

Proof:  We  first  show  that  if  u  is  dequeued  before  v,  then 
D[u]  <  D [t;]  at  the  instants  when  they  are  dequeued.  We  prove 
this  by  induction.  The  initial  condition  holds  trivially.  Then,  as¬ 
sume  it  is  true  for  the  first  I  dequeued  vertices  zi, ...  ,zi.  Con¬ 
sider  the  (/  +  l)th  dequeued  vertex  .  At  the  instant  when  zi  is 
dequeued,  if  D[zi^i]  is  updated  based  on  D[zi]  in  Step  2,  then 
D[zi]  <  D[zi^i]  even  after  the  update.  If,  on  the  other  hand, 
D[zi^i]  is  not  updated  at  this  instant,  then  D[zi]  <  D[zi^i] 
given  that  the  dequeued  vertex  has  the  smallest  distance. 

Next,  we  prove  the  lemma  by  induction.  From  Step  1  of  the 
algorithm,  all  the  affected  vertices  v  in  e  will  be  dequeued  first 
with  E[v]  =  6,  P[v]  =  X,  and  D[v]  =  d'[x]  w{e).  Based  on 

Lemma  2,  D[v]  <  f[u]  +  w{e)  for  any  u  e  e.  It  thus  follows 
that  the  hyperpath  to  v  through  x  and  e  is  the  shortest  one  with 
D[v]  =  d'[v]. 

Assume  for  zi^ ...  ^zi,  D[zi]  =  f[zi]  are  satisfied  for  all 
i  =  1, . . . ,  /.  Consider  the  (/  +  l)th  dequeued  vertex  zi^i  ^  e. 
Let  u  =  P[zi^i]  be  its  parent  in  the  new  shortest  hyperpath. 
Then,  based  on  the  fact  that  distances  of  the  dequeued  vertices 
are  monotonically  increasing  with  the  order  of  the  dequeueing 
as  shown  at  the  beginning  of  the  proof,  u  cannot  be  any  vertex 
dequeued  after  zi^i.  Since  zi^i  0  e,  it  is  also  clear  that  u  cannot 
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be  an  unaffected  vertex  (otherwise,  zi^i  will  be  unaffected, 
which  contradicts  Lemma  3).  We  thus  have  u  G  . . .  ,zi}. 
Let  u  =  Then,  when  Zi  is  dequeued,  D[zi^i]  will  be  updated 
to  the  shortest  distance  d'[zi^i]  due  to  the  induction  hypothesis 
of  D[zi]  =  d'[zi].  This  completes  the  proof.  ■ 

Based  on  Lemmas  3  and  4,  the  shortest  distances  of  all  af¬ 
fected  vertices  will  be  updated  correctly.  Based  on  Lemma  3, 
no  unaffected  vertex  will  be  enqueued,  and  their  distances  re¬ 
main  the  same.  It  is  not  difficult  to  see  from  the  algorithm  that 
P[v]  and  E[v]  are  also  correctly  maintained  for  all  v. 


Appendix  B 
Proof  of  Theorem  2 

We  first  show  the  correctness  of  the  coloring  process  as  given 
in  the  following  lemma. 

Lemma  5:  The  coloring  process  correctly  colors  all  the  af¬ 
fected  vertices. 

Proof:  We  first  state  the  following  simple  facts  without 
proof:  Given  a  relationship  tree,  after  the  hyperedge  weight  in¬ 
crease:  1)  if  TJ  is  pink  or  white,  then  all  its  descendent  in  this 
relationship  tree  are  white;  2)  if  a  ?;  is  red,  then  all  its  children 
in  the  relationship  tree  are  either  pink  or  red;  3)  if  u  is  affected, 
either  v  £  e  ox  P[v\  is  red.  These  facts  can  be  directly  obtained 
from  the  definition  of  the  colors.  It  is  also  easy  to  see  that  ver¬ 
tices  are  dequeued  from  M  in  a  nondecreasing  order  of  their 
current  distance  D[-] .  This  is  because  each  time  a  vertex  2;  is  de¬ 
queued  from  M,  the  possible  new  vertices  to  be  enqueued  into 
M  are  2;  ’s  children  with  distances  no  smaller  than  D  [2;] . 

Then,  the  proof  of  the  lemma  has  two  parts:  First,  we  prove 
that  all  affected  vertices  are  enqueued  into  M ;  then  we  prove  by 
induction  that  only  affected  vertices  are  enqueued  into  M  and 
their  colors  are  correctly  identified. 

We  prove  the  first  part  by  contradiction.  Assume  that  there 
exists  an  affected  vertex  v  that  is  not  enqueued  into  M.  It  is 
easy  to  see  that  v  ^  e  because  all  the  affected  vertices  in  e  are 
enqueued  in  Step  1 .  Based  on  the  third  fact  stated  above,  P  [tj]  is 
red.  Based  on  the  hypothesis,  P[v]  is  not  enqueued  (otherwise,  v 
will  be  enqueued  in  Step  2).  Continuing  this  line  of  arguments, 
we  eventually  reach  the  root  of  the  relationship  tree  and  arrive 
at  the  contradiction  that  the  source  s  is  red. 

We  prove  the  second  part  by  induction.  It  is  easy  to  see  that 
all  the  vertices  initially  enqueued  into  M  are  affected  vertices. 
It  remains  to  show  that  the  first  vertex  2^1  dequeued  from  M  is 
colored  (pink  or  red)  correctly.  To  show  that,  we  need  to  estab¬ 
lish  that  the  algorithm  correctly  determines  whether  there  is  an 
alternative  shortest  hyperpath  to  2^1  with  the  same  distance,  i.e., 
d[zi]  =  f[zi].  The  key  here  is  to  show  that  checking  the  cur¬ 
rently  nonred  neighbors  (which  may  become  red  in  the  future) 
of  2^1  will  not  lead  to  a  false  alternative  path.  This  follows  from 
the  fact  that  zi  has  the  smallest  distance  D[‘]  among  all  affected 
vertices  (which  belong  to  the  set  of  vertices  consisting  of  the 
affected  vertices  in  e  and  their  descendents). 

Next,  assume  that  vertices  zi,  Z2, . . . ,  zi  dequeued  from  M 
are  all  affected  vertices  and  are  correctly  colored.  Consider  the 
next  dequeued  vertex  .  It  is  an  affected  vertex  because  it  is 


either  enqueued  in  Step  1  with  E[v]  —  e  or  enqueued  in  Step  2 
with  a  red  parent.  To  show  that  zij^i  will  be  colored  correctly, 
we  use  a  similar  argument  by  showing  that  the  currently  nonred 
neighbors  of  zij^i  will  not  give  a  false  alternative  path.  The  latter 
follows  from  the  fact  that  all  affected  vertices  will  be  enqueued 
and  those  dequeued  after  zi^i  have  distances  no  smaller  than 
D[zi^i].  This  completes  the  induction.  ■ 

We  now  show  that  D[v]j  P[v]  and  E[v]  are  correctly  main¬ 
tained  for  all  u.  For  each  red  vertex  v,  its  distance  is  set  based  on 
the  current  shortest  distance  from  a  nonred  neighbor  in  Step  3. a. 
The  rest  of  the  algorithm  is  essentially  Gallo’s  extension  of  Di- 
jkastra’s  algorithm  with  the  current  initial  distance.  The  cor¬ 
rectness  of  the  algorithm  thus  follows.  It  is  not  difficult  to  see 
that  P[-]  and  E[']  are  correctly  updated  for  both  red  and  pink 
vertices. 
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